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ABSTRACT 

Quasar host galaxies are key for understanding the relation between galaxies and the supermassive 
black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host 
galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky 
Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable 
decomposition of the observed spectra into nuclear and host components, using spectral models of 
quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate 
age, mass M*, and velocity dispersion tr* of the host stars, the star formation rate (SER), quasar 
luminosity, and SMBH mass M,, for each object. The quasars are preferentially hosted by massive 
galaxies with M* ~ 10®characterized by stellar ages around a billion years, which coincides with 
the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies 
have relatively low SERs and fall below the main sequence of star-forming galaxies at similar redshifts. 
These facts suggest that the hosts have experienced an episode of major star formation sometime in 
the past billion years, which was subsequently quenched or suppressed. The derived M, — cr* and 
M, — M* relations agree with our past measurements and are consistent with no evolution from the 
local Universe. The present analysis demonstrates that reliable measurements of stellar properties of 
quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers 
with future powerful instruments such as the Subaru Prime Eocus Spectrograph. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: nuclei — galaxies: stellar content 
— quasars: general — quasars: supermassive black holes 


1. INTRODUCTION 

Over the last few decades, it has become clear that su¬ 
permassive black holes (SMBHs) are ubiquitous at galaxy 
centers throughout the Universe. Stellar or gaseous 
motions at galactic centers indicate that almost every 
galaxy bulge in the local Universe harbors a SMBH 
(e.g.. Dressier 1989; Kormendy 1993; Magorrian et al. 
1998). The correlation between SMBH mass (A/,) and 
the stellar velocity dispersion (cr*) or mass of the bulges 
(e.g., Merritt & Ferrarese 2001; McLure & Dunlop 2002; 
Haring & Rix 2004; Giiltekin et al. 2009) suggests that 
the two coevolve, or at least strongly influence each other. 
The observed scatter in the M, — cr* relation is remark¬ 
ably small (Ferrarese & Merritt 2000; Gebhardt et al. 
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2000). This situation may require a fine-tuning pro¬ 
cess, which is perhaps a combination of internal physics 
within galaxies or dark halos and mass averaging via fre¬ 
quent mergers (see the recent review by Kormendy & Ho 
2013, and references therein). Beyond the local Universe, 
SMBHs are known to exist even at very high redshifts. 
Recent wide-held surveys have discovered more than 40 
quasars at z > 6 (e.g.. Fan et al. 2006; Jiang et al. 2009; 
Willott et al. 2010; Bafiados et al. 2014), with the most 
distant discovery at z = 7.1 (Mortlock et al. 2011). The 
estimated SMBH mass exceeds M, = 10®Mq in many 
cases (e.g., Jiang et al. 2007; Wu et al. 2015), which has 
initiated intense discussion on the process of their forma¬ 
tion and evolution in the short time scales of the early 
Universe (e.g., Ferrara et al. 2014; Madau et al. 2014). 
The M, — tr, relation of these objects tend to be offset 
toward higher M, compared to the local relation at hxed 
cr* (e.g., Wang et al. 2010), but possible biases in sample 
selection and measurements complicate the interpreta¬ 
tion (see, e.g., Schulze & Wisotzki 2011; Shen 2013). 

Since the Soltan (1982) argument and the latest obser¬ 
vations suggest that the local SMBHs have acquired most 
of their mass during an active galactic nucleus (AGN) 
phase (e.g., Yu & Tremaine 2002; Marconi et al. 2004; 
Brandt & Alexander 2015; Gomastri et al. 2015), a key 
to understanding SMBH growth and the origin of the 
M, — (7, relation may be found in observational proper¬ 
ties of AGNs. One of the outstanding issues is the nature 
of their host galaxies; in which galaxies do AGNs occur, 
and what impact do they have on the galaxies? These 
questions have been addressed by numerous studies in re- 
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cent years, especially in terms of negative feedback effects 
of AGNs on the star formation rates (SFRs) of the host 
galaxies. It has been suggested that the energy and/or 
momentum input of AGNs, which may be triggered along 
with starbursts by interactions/mergers of gas-rich galax¬ 
ies (e.g., Barnes & Hernquist 1991), could suppress star 
formation by expelling the cold gas from galaxies (e.g., Di 
Matteo et al. 2005; Springel et al. 2005) and/or heating 
the gas in dark halos (e.g., McNamara & Nulsen 2007). 
Indeed, recent observations of AGN-driven gas outflows 
on small scales (e.g.. Pounds et al. 2003; Trump et al. 
2006; Ganguly et al. 2007; Reeves et al. 2009; Tombesi et 
al. 2010; Gofford et al. 2013) as well as large scales (e.g., 
Nesvadba et al. 2006; Feruglio et al. 2010; Greene et al. 
2011; Cano-Diaz et al. 2012; Liu et al. 2013a,b) and ex¬ 
tended emission-line regions photoionized by AGNs (e.g., 
Fu & Stockton 2009; Husemann et al. 2010; Matsuoka 
2012; Keel et al. 2012) strongly suggest that AGNs have 
a significant impact on their host galaxies. The idea of 
negative AGN feedback is also favored because it may 
provide a solution to longstanding problems with galaxy 
formation models; for example, it is currently the most 
compelling process to reconcile the quite different shapes 
of the dark halo mass function predicted in the AGDM 
cosmology and the observed galaxy stellar mass function 
at the high mass end (e.g., Bower et al. 2006; Groton et 
al. 2006; Somerville et al. 2008; Dave et al. 2011; Ghoi 
& Nagamine 2011). The AGN feedback may have two 
distinct modes, i.e., quasar mode and radio mode, which 
become efficient at different phases of SMBH growth and 
galaxy evolution (see, e.g., Fabian 2012). 

There have been a number of studies on host galax¬ 
ies of various AGN classes. Different types of AGNs 
may play different roles in galaxy evolution (e.g., Hickox 
et al. 2009; DiPompeo et al. 2014, 2015), hence it is 
crucial to understand the entire AGN and host popula¬ 
tions. Elitzur (2012) pointed out that optically-classified 
narrow-line AGNs are more likely found in dusty cir- 
cumnuclear environments than are broad-line AGNs (see 
also Netzer 2015). In a merger-driven scenario of galaxy 
and SMBH co-evolution (Hopkins et al. 2006), a major 
merger Hrst creates the dusty starburst/AGN phase ob¬ 
served as luminous infrared galaxies, and then evolves 
into the dust-free phase observed as optically-luminous 
quasars (see also Sanders et al. 1988). Therefore, AGNs 
selected at different wavelengths may have fundamentally 
different host properties; for instance, dusty infrared- 
selected AGNs (e.g., Lacy et al. 2013) may be associ¬ 
ated with more active star formation on average than 
are optically-selected AGNs. 

Early studies on this subject (e.g., Boroson & Oke 
1982; Kotilainen & Ward 1994; Ronnback et al. 1996; 
Bahcall et al. 1997; McLure et al. 1999; Kirhakos et al. 
1999; Dunlop et al. 2003; Jahnke et al. 2004, 2007) have 
shown that host galaxies of optically-luminous AGNs are 
not a random subsample of normal galaxies. Instead, 
these AGNs were found to reside in luminous galaxies 
with young stars and a higher than usual fraction of close 
companions or tidal features. It has also been suggested 
that the host galaxies have a substantial contribution 
from intermediate-age stars (e.g., Canalize & Stockton 
2013), which may correspond to the late stage of accre¬ 
tion episode predicted in the above merger-driven evolu¬ 
tionary scenario. Recent studies have taken advantage of 


a large sample of obscured AGNs, which show no visible 
sign of accretion-disk emission or broad lines at ultra¬ 
violet (UV) to optical wavelengths, discovered in large 
spectroscopic surveys (e.g., Kauffmann et al. 2003; Za- 
kamska et al. 2003; Hao et al. 2005; Reyes et al. 2008) or 
X-ray observations (e.g., Nandra et al. 2007; Georgakakis 
et al. 2008; Silverman et al. 2008; Mullaney et al. 2012). 
Dust extinction serves as a natural coronagraph in these 
objects and makes analysis of the host galaxies much 
easier than in unobscured populations. Overall, these 
studies have statistically confirmed the earlier findings 
and shown that obscured AGNs reside preferentially in 
massive galaxies. While there have been controversies 
over the host colors, recent work points out the impor¬ 
tance of using a mass-matched control sample of inactive 
galaxies for unbiased measurements. Once the mass de¬ 
pendence is properly accounted for, the AGN fraction is 
more or less constant across the host colors or moder¬ 
ately enhanced in bluer galaxies (e.g., Silverman et al. 
2009; Xue et al. 2010; Hainline et al. 2012; Rosario et al. 
2013; Trump et al. 2015). At the same time, it has been 
demonstrated that the color by itself may not be a good 
tracer of the star-forming properties of the host galaxies 
(e.g., Gardamone et al. 2010; Rosario et al. 2013). 

As a complementary effort to the studies of obscured 
AGNs in large contemporary surveys, we are investigat¬ 
ing the host galaxies of “classical” broad-line quasars 
selected at optical wavelengths. In our previous work 
(Matsuoka et al. 2014a), we analyzed the images of about 
1,000 Sloan Digital Sky Survey (SDSS; York et al. 2000; 
Eisenstein et al. 2011) broad-line quasars at 2 : < 0.6 and 
explored their host properties. Deep co-added SDSS im¬ 
ages of the Stripe 82 region (Annis et al. 2014) in five 
bands {u, g, r, i, and z; Eukugita et al. 1996) were 
used for the analysis. We developed a technique to de¬ 
compose spatially a quasar image into nuclear and host 
components, using the point spread function (PSE) and 
Sersic (1968) models, and measured the colors and stellar 
masses of the host galaxies. We showed that the quasars 
are hosted exclusively by massive galaxies and that these 
host galaxies are considerably bluer than the red se¬ 
quence, showing stark contrast to the color-magnitude 
diagram (CMD) of inactive galaxies (see also Trump et 
al. 2013). We also found a positive correlation between 
M, and stellar mass M*, but the relation is offset to¬ 
ward large M, compared to the local relation at fixed 
M*. While this result could indicate that SMBHs grow 
earlier than do the host galaxies, we deferred a definitive 
conclusion until observational biases are fully understood 
(e.g., Schulze & Wisotzki 2011; Shen 2013). 

In this paper, we present the results of a spectral 
decomposition analysis of SDSS broad-line quasars at 
2 ; < 1. By making use of the high signal-to-noise ra¬ 
tio (SNR) spectra produced by the SDSS Reverberation 
Mapping (RM) project (see the following section), we 
measure the spectral properties of about 200 quasar host 
galaxies and explore their nature. This paper is orga¬ 
nized as follows. The data and sample are presented 
in §2. §3 describes the technique used to decompose a 
quasar spectrum into nuclear and host components, and 
the associated error assessments. The main results ap¬ 
pear in §4 and are discussed in context of galaxy and 
SMBH evolution in §5. The summary follows in §6. 
The cosmological parameters Hq = 70 km s“^ Mpc“^, 
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Hm = 0.3, and Ha = 0.7 are used throughout this work. 
All magnitudes are presented on the AB system (Oke & 
Gunn 1983). 

2. DATA ACQUISITION AND BASIC PROCESSING 

We use the data acquired in the SDSS-RM project, 
in which a single spectroscopic field was repeatedly ob¬ 
served to explore the variability of quasars. The full 
technical details of the project are found in Shen et 
al. (2015a), and hence we repeat only the most rel¬ 
evant aspects here. SDSS-RM was conducted during 
the dark/grey time in the final season (2013-2014) of 
the SDSS-III Baryon Oscillation Spectroscopic Survey 
(BOSS; Dawson et al. 2013) following its early com¬ 
pletion due to unexpectedly good weather. The tar¬ 
get field (centered at aj 2 ooo = 14*'14“49®.00, <5j2ooo 
= -|-53°05™00®.0) coincides with the Panoramic Survey 
Telescope & Rapid Response System 1 (Pan-STARRS 1; 
Kaiser et al. 2010) Medium Deep Field MD07 which lies 
within the Canada-France-Hawaii Telescope Legacy Sur¬ 
vey W3 field. About 1,000 spectroscopically confirmed 
quasars are known in this 7-deg^ field from SDSS-I/II 
and SDSS-III, most of which were targeted by the BOSS 
spectroscopy. Additional quasars have been identified 
by the DEEP2 Galaxy Redshift Survey (Newman et al. 
2013) and by follow-up spectroscopy of Pan-STARRS 
1 variable sources (P. J. Green et al., in preparation). 
The redshifts of these quasars have been determined by 
matching identified emission lines with a pre-defined line 
list or by cross correlating observed spectra with spectral 
templates. 

The SDSS-RM quasars were selected as a flux-limited 
sample of 849 objects at i < 21.7 mag, with objects in 
fiber collisions being removed. Each quasar is labelled 
with an identification number RMID. We stress that the 
entire sample is composed of spectroscopically-confirmed 
broad-line quasars, the vast majority of which were se¬ 
lected at optical wavelengths; further details of the sam¬ 
ple selection are given in §3.1 of Shen et al. (2015a). 
The quasars follow the luminosity-redshift distribution 
expected from the luminosity function estimates in Hop¬ 
kins et al. (2007), suggesting that the sample is fairly 
uniformly selected above the adopted flux limit. Figure 
1 displays the distribution of their redshifts and abso¬ 
lute magnitudes in SDSS i band (Mi). Mi was derived 
from the SDSS PSF magnitudes, /c-corrected to z = 0 
with the correction factors taken from Richards et al. 
(2006). The quasars analyzed in this work (z < 1; see 
§3.2) have luminosities comparable to or lower than the 
classical quasar threshold Mb = —23 mag (Schmidt & 
Green 1983), which corresponds to Mi ~ —23 mag as¬ 
suming the SDSS quasar composite spectrum of Vanden 
Berk et al. (2001) and correcting for the different cosmol¬ 
ogy used in Schmidt & Green (1983). 

The SDSS-RM observations were carried out with the 
BOSS spectrograph (Smee et al. 2013) mounted on the 
SDSS 2.5-m telescope (Gunn et al. 2006) in 32 epochs 
from January to July 2014. The spectrograph provides 
wavelength coverage from A = 3,650 to 10,400 A with 
a spectral resolution of i? 2,000. The BOSS plates 
allow 1,000 fibers of 2" diameter to be observed simulta¬ 
neously; the fibers were allocated to the 849 quasars, 70 
standard stars, a luminous red galaxy (required by the 
BOSS targeting algorithm), and 80 sky positions. A typ- 



Fig. 1. — Redshifts and z-band absolute magnitudes of the 849 
SDSS-RM quasars (small and large dots), defined by the flux limit 
i < 21.7 mag (dashed line). We analyze 191 quasars at z < 1.0 
(dotted line) in this study. Our final sample consists of 156 quasars 
(large dots) with successful spectral decomposition (see §3.2). 

ical epoch consists of eight to ten 15-minute exposures, 
totaling 65 hours when summed over the 32 epochs. The 
data were first processed through the BOSS spectro¬ 
scopic pipeline idlspec2d v5_7_l, which performs flat- 
fielding, Id extraction, wavelength calibration, sky sub¬ 
traction, and flux calibration (Bolton et al. 201^ and 
references therein). The spectrum was resampled to a 
pixel scale of 69 km s“^, which corresponds to roughly 
a half of the resolution element. The flux calibration 
is tied to the PSF magnitudes of standard stars in the 
plate. These pipeline-processed data were released to 
the public in the SDSS data release (DR) 12 (Alam et 
al. 2015). The spectrophotometry was further improved 
with a custom flux-calibration algorithm; SDSS-RM ob¬ 
served 3.5 times more standard stars per plate than the 
main BOSS survey, which enabled us to trace the spatial 
variation of the calibration more precisely (Shen et al. 
2015a). Later in this work, we apply an additional aper¬ 
ture correction so that the decomposed host magnitude 
represents the total galaxy light. The data of all epochs 
were stacked by inverse-variance weighted mean to cre¬ 
ate high SNR spectra, which are used throughout this 
work. The SNR of the stacked spectra range from ~5 to 
^^300 per pixel, with the median value of ^30. The pixel 
flux errors were taken from the BOSS pipeline and prop¬ 
erly propagated throughout the post-processing. We cor¬ 
rected the spectra for Galactic dust extinction with the 
extinction curve of Pei (1992) and the color excess E-q-y 
taken from Schlegel et al. (1998). 

3. DATA ANALYSIS 

We describe the details of the data analysis in this sec¬ 
tion. In §3.1, we introduce our spectral decomposition 
technique which separates the galaxy component from 
an observed quasar spectrum. This method is applied to 
the SDSS-RM data in §3.2. In §3.3, we evaluate the relia¬ 
bility of the present method by Monte Carlo simulations 
and by comparisons with previous measurements. 
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3.1. Spectral decomposition technique 

We aim to extract host-galaxy information from the 
observed spectra with as little ambiguity as possible. 
This is a challenge, since different combinations of quasar 
and stellar spectral models give rise to similar total spec¬ 
tra at UV-to-optical wavelengths. For the present anal¬ 
ysis, we choose the rest-frame spectral window from A 
= 3700 A to 5400 A, which includes a wealth of stellar 
absorption features as well as important nebular lines 
such as [O II] A3727, H/3, and [O III] AA4959, 5007. 
We devise a method to apply a series of spectral fits 
to the data in such a way that degeneracies between 
model components are minimized. The following models 
are used to represent the relevant emission components, 
namely, the quasar accretion disk, gas in the broad line 
region (BLR), narrow line region (NLR), and interstellar 
medium (ISM), and stars in the host galaxy. 

1. Accretion disk: the emission from the quasar ac¬ 
cretion disk can be represented by a single power- 
law continuum {f\ oc A“p') in our spectral window, 
with the amplitude and slope being free parame¬ 
ters. Previous observations have found slopes span¬ 
ning the range of —2.5 < Opi <0.0, with a typical 
value of Opi ~ —1.5 (Richstone & Schmidt 1980; 
Francis 1996; Vanden Berk et al. 2001; Ivezic et al. 
2002; Pentericci et al. 2003). 

2. Gas in the BLR, NLR, and ISM: 

(a) 3700 A < A < 4450 A: numerous high- 
order Balmer lines and metal lines as well 
as the Balmer continuum make up a pseudo¬ 
continuum in this spectral region. It is mod¬ 
eled by subtracting the power-law contribu¬ 
tion from the SDSS quasar composite spec¬ 
trum provided by Shen et al. (2011). We take 
the composite of the highest-luminosity bin 
(Lsioo = 1045-5 erg s"!) so that the host con¬ 
tamination is minimized; Shen et al. (2011) 
report that the contamination at A = 5100 A 
is negligible in a quasar with L 5100 > 1045-5 
erg s“5. The power-law contribution is deter¬ 
mined in two continuum windows (in which 
strong emission lines are absent) at A = 4200 
- 4230 and 5470 - 5500 A. This pseudo¬ 
continuum template also includes individual 
BLR, NLR, and ISM lines present in the orig¬ 
inal composite spectrum, although the ISM 
contribution should be minimal at this highest 
quasar luminosity. We consider seven addi¬ 
tional narrow lines or line complexes, namely, 
[O II] A3727, [Fe VII] A3760, [Ne III] A3870, 
He I -f H 8 A3890, [Ne III] A3969, Hd, and H 7 , 
to fit the ISM lines. They are each modeled 
with Gaussian profiles sharing the same veloc¬ 
ity dispersion (ag) and velocity offset relative 
to the object’s formal redshift (u°®), both of 
which are free parameters. 

(b) 4450 A < A < 5400 A: the dominant con¬ 
tributors to this spectral region are the Fe II 
pseudo-continuum, H/3, and [O III] A4959, 
A5007 lines. We use the empirical Fe II 


template of Veron-Cetty et al. (2004) cre¬ 
ated from an observed spectrum of I Zw 1, a 
prototype narrow-line Seyfert 1 galaxy (e.g., 
Osterbrock & Pogge 1985). This template 
is velocity-broadened with a Gaussian kernel 
with CTpeii = 1,500 km s“4; we test a differ¬ 
ent value of CTpeii later. H/3 is modeled with 
three Gaussian emission profiles with variable 
amplitudes, widths, and velocity offsets, while 
the two [O III] lines are each modeled with 
two Gaussians. The widths of these lines 
are varied independently from those of the 
above ISM lines. We assume that the two 
[O III] lines have the theoretical intensity ra¬ 
tio A5007/A4959 = 3.0 and that the two lines 
share the same profile (i.e., widths and veloc¬ 
ity offsets). 

3. Stars: we use simple stellar population (SSP) mod¬ 
els of Maraston & Strbmback (2011) to represent 
the stellar emission. SSP is a single generation of 
stars formed in an instantaneous starburst; it is the 
simplest assumption for a stellar population and is 
suitable for the present analysis, in which the dom¬ 
inant quasar light makes it difficult to exploit full 
details of the stellar properties. We test more com¬ 
plicated stellar population models in §3.3. For SSP 
models, we adopt the STELIB (Le Borgne et al. 
2003) stellar spectral library, the Ghabrier (2003) 
initial mass function (IMF), and solar metallicity. 
Stellar age (G) takes values from 30 Myr to the age 
of the Universe at the redshift of each object. The 
SSP spectra are velocity-broadened in accordance 
with the stellar velocity dispersion ) and the in¬ 
strumental resolution^5 assumed to be (Tinst = 65 
km s“5^ and offset in wavelength by relative 
to the object’s redshift; both cr* and are free 
parameters. The amplitude of the SSP model de¬ 
termines the stellar mass M*. We test different 
choices of stellar library, IMF, and metallicity in 

§4- 

Our default models do not incorporate the effect of 
dust extinction. In fact, the amount of extinction ob¬ 
served in optically-selected broad-line quasars is usually 
quite small (Richards et al. 2003; Hopkins et al. 2004; 
Salvato et al. 2009; Matute et al. 2012) and its effect is 
limited. Recently Krawczyk et al. (2014) investigated a 
large sample of SDSS broad-line quasars and found that 
the vast majority are consistent with Eb-v < 0.1 mag. 
Here we assume no dust extinction for both quasar nu¬ 
clei and host galaxies, and later conhrm that non-zero 
Eb-v values have little impact on our qualitative con¬ 
clusions. We note that there are dust-reddened broad¬ 
line quasars which can be identified in the near-infrared, 
although their surface density is much lower than that 
of optically-selected quasars (e.g., Glikman et al. 2012). 
Urrutia et al. (2008) presented an enhanced merger frac¬ 
tion in dust-reddened quasars, hence their host stellar 
properties may be quite different from those of normal 
optically-selected quasars. 

The actual instrumental resolution varies as a function of 
wavelength, but here we assume the fixed mean value (see Smee et 
al. 2013). 
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Fig. 2. — Examples of the spectral decomposition. Each row corresponds to a single quasar in the full spectral window (left) and expanded 
view (right). The first, second, and third rows show objects fitted with young (t* = 30 Myr), intermediate (1 Gyr), and old (8 Gyr) SSP 
models, respectively, while the last row displays an object with strong narrow ISM emission lines. Black lines indicate the observed co-added 
spectra, while the red lines show the total model spectra. The other lines represent the individual model components, namely, the quasar 
power-law continuum (blue dashed), the gas emission from the BLR, NLR, and ISM (blue solid), and the stellar emission (green). 


We perform a series of spectral fits rather than fitting 
all the model components simultaneously. First, we de¬ 
termine the stellar age t* by fitting the spectrum over 
the rest-frame wavelength range A = 3700 - 4050 A with 
the quasar power-law slope fixed to Opi = —1.5. This 
spectral window is chosen because (i) a wealth of strong 
stellar absorption features are present, including the Ca 
HK lines and the 4000 A break; (ii) the quasar emission 
is relatively featureless; (iii) the precise value of api is not 
important, since the flux ratio between the two edges of 
the window (/ 4050 // 3700 ) changes only by ~10 % when 
api varies from —1.5 to the extreme values of 0.0 or —2.5. 
Second, the wider spectral range of A =3700 - 4750 A 


and 5050 - 5400 A (i.e., the full 3700 - 5400 A spectral 
window excluding the H/3 region) is fit by varying all pa¬ 
rameters but t*, which is fixed to the value determined 
above. We then return to the first step and update t* 
with the latest app These two steps are iterated until t* 
converges, but in most cases the convergence occurs with 
no more than two iterations. Finally, we subtract the 
best-fit model from the full spectrum and fit the residual 
with the H/? and [O III] line models in the range A = 
4750 - 5050 A. 
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3.2. Application to the data 

We apply the above algorithm to all the 191 SDSS- 
RM quasars found at 2 < 1; our spectral window (A = 
3700 - 5400 A) is shifted beyond the spectrograph cov¬ 
erage at higher redshifts. The standard method is 
used for the model fitting with the IDL routine MPFIT 
(Markwardt 2009). As described above, the pixel flux 
errors were taken from the BOSS spectroscopic pipeline 
and propagated appropriately through the co-addition 
process. Figure 2 shows four typical examples of the de¬ 
composition results. While the spectra are usually domi¬ 
nated by quasar light, most of the small-scale features are 
fit by the stellar models. The median reduced of the 
fits over the whole spectral window is about five; as seen 
in Figure 2, this rather large value is caused by the sys¬ 
tematic discrepancy at certain wavelengths between the 
observed spectra and the models, the latter being not 
complex enough to provide perfect fits to every emission 
and absorption feature across the spectrum. The uncer¬ 
tainties in the best-fit parameters are estimated following 
the Monte Carlo method outlined in Shen et al. (2011). 
We create 50 mock spectra by adding random noise to the 
original spectrum, using the Gaussian probability density 
function (PDF) with the standard deviation equal to the 
flux error at each pixel, and process them through the 
decomposition algorithm. The measurement errors are 
evaluated from the 68 % central intervals of the param¬ 
eter distributions of the 50 trials. We add additional 20 
% error to the derived stellar age, which accounts for the 
minimum uncertainty arising from non-continuous cov¬ 
erage of ages in the SSP models. We perform further as¬ 
sessment of possible systematic uncertainties in the next 
section. 

Figure 3 presents the histogram of the fractional con¬ 
tribution of the stellar component to the total flux (/*) 
measured at A = 4000 A in the rest frame. We exclude 
the 31 objects with /* < 0.1 from the following anal¬ 
ysis; we will demonstrate in the next section that their 
host results are unreliable. An additional four objects are 
excluded because they have particularly large statistical 
uncertainties in the host stellar mass, A\og{M^:/M q) > 
0.5 dex. These excluded objects tend to have high quasar 
luminosities, as shown in Figure 1. Our final sample con¬ 
sists of the remaining 156 quasars, more than 80 % of 
the initial 191 objects. The median redshift and abso¬ 
lute magnitude of the 156 quasars are {z) = 0.72 and 
(Mi) = —22.2 mag, respectively. 

We focus on six physical quantities in this work, 
namely, stellar age O, stellar mass M*, SFR, stellar ve¬ 
locity dispersion cr*, [O III] A5007 luminosity T[oiiip and 
SMBH mass M,. The parameters O, M*, a*, and A[oiii] 
are taken directly from the decomposition results. Since 
the quantity ct* is poorly constrained in many cases, its 
measurement is rejected if the best-fit value lies outside 
the physically plausible range (10 < cr* < 400 km s“^) or 
smaller than three times the estimated statistical error. 
SFRs are estimated from the [O II] A3727 luminosity. 
As shown in Figure 2, the [O II] line is comprised of two 
components, i.e., the quasar BLR/NLR template and the 
host ISM line with Gaussian profile. The former contri¬ 
bution is tied to the overall amplitude of the BLR/NLR 
template, which is fit over the full spectral window. We 
use the luminosity of the latter ISM component (Ljoiij) 



Stellor fraction (fjt^r) 


Fig. 3.— Histogram of the fractional stellar contribution to the 
total flux (/*), measured at A = 4000 A in the rest frame, for 
the 191 SDSS-RM quasars at .2 < 1. The typical /* uncertainty 
is shown with the error bar at the top right corner. We deem 
the decomposition analysis successful for 156 objects (82 % of the 
initial 191 quasars) with /* > 0.1 (dotted line). 


for SFR estimates, assuming the Kennicutt (1998) cali¬ 
bration: 

SFR (M© yr-i) = (1.4 ±0.4) x IQ-^^ Ljo„] (erg s”'). 

( 1 ) 

Since the NLR usually contributes only weakly to this 
low-ionization line, it may be a good tracer of star for¬ 
mation in quasar host galaxies even without spectral de¬ 
composition (Ho 2005). On the other hand, Kim et al. 
(2006) found that the [O II]/[O III] ratios observed in 
SDSS AGNs at 2 < 0.3 are fully consistent with AGN 
photoionization, hence there is no need to invoke any ad¬ 
ditional [O II] source such as star formation. Although 
the mean quasar contribution is properly subtracted with 
the BLR/NLR template in this work, we conservatively 
treat the derived SFRs as rough estimates only, given 
that our single BLR/NLR template cannot trace object- 
to-object variation of the [O II] strength. Another con¬ 
cern is the effect of dust extinction, which we discuss 
later. 

SMBH masses are estimated with the single-epoch 
virial estimator (Vestergaard & Peterson 2006): 


log 




FWHM(H/3) 


1000 km s“i 
±(6.91 ±0.02). 


A£a(5100 A) 
lO"^'^ erg s“i 


0.5\ 


( 2 ) 


It is believed that this H/3-based estimator is more re¬ 
liable than other single-epoch estimators, such as those 
based on G IV A1549 or Mg II A 2800 (Shen & Liu 2012; 
Shen 2013). In some cases our H/1 models apparently fit 
the underlying continuum residuals rather than the line 
itself, especially when the SNR around H/3 is low. The 
M, measurements were rejected for 16 such fits identi¬ 
fied by visual inspection. We will eventually obtain more 
precise SMBH masses for all the objects with RM anal¬ 
ysis. 

The limited aperture size of the BOSS spectrograph 






Host galaxies of SDSS-RM quasars 


7 



40.0 40.5 41.0 41.5 42.0 42.5 43.0 7.6 7.8 8.0 8.2 8.4 8.6 8.8 9.0 

logCLfoiiil/erg s"') loglMsH/Msun) 


Fig. 4. — Deviations of the measured (output) values from the intrinsic (input; the quantities with the superscript “0”) values for (top 
left), M* (top right), SFR (middle left), tr* (middle right), F[oiii] (bottom left), and M, (bottom right) for simulated spectra. The black 
dots represent the spectra with reliable fits of the stellar component (output stellar fraction /* > 0.1, mass error Alog M, < 0.5 dex, and 
stellar age t, > 0.5 Gyr), while the blue dots represent those with younger stellar ages (/« > 0.1, Alog M« < 0.5 dex, and t* < 0.5 Gyr; we 
find few such cases in the real quasars as described in the text). The open circles show the remaining spectra with poor fits of the stellar 
component (/, < 0.1 or Alog M« > 0.5 dex). The symbols in the upper left panel are given small random offsets in stellar age to improve 
visibility. The typical statistical error is shown by the error bar at the bottom left corner of each panel. The horizontal lines represent the 
identity lines. The numbers in each panel indicate, for the black dots, the RMS scatter around the identity line and the fraction (/id) of 
those whose ordinate values are consistent with zero within Icr statistical error. 


fibers does not include all the light from the extended 
hosts. The fiber diameter of 2" corresponds to approx¬ 
imately 14 kpc at the median redshift of our sample 
(z = 0.72). We correct the stellar mass, derived from 
the decomposed galaxy flux (/g*’®’’), for the aperture loss 
as follows. Since the spectrophotometry of the BOSS 
spectroscopic pipeline is tied to the PSF magnitude of 
standard stars, the quasar flux decomposed from each 
spectrum (/q) should correctly represent the total flux 
of the unresolved quasar nucleus. The SDSS imaging 
pipeline measures the total nuclear plus host flux (/tot) 
with the cModel algorithm (Abazajian et al. 2004), which 
fits galaxy models to the observed image profile. There¬ 


fore, the total host flux (/g) is given by fg = /tot — /q- 
Multiplying the stellar mass of each source measured 
from the spectra by the ratio fg/fg^^'^ provides a crude 
estimate of the total stellar mass, with the implicit as¬ 
sumptions that the mass-to-luminosity ratio in the fiber 
represents that in the whole galaxy and that quasar vari¬ 
ability between the epochs of imaging (before 2005) and 
spectroscopic (2014) observations does not introduce a 
significant bias. This aperture correction is calculated 
in the observed-frame r band. The correction is not ap¬ 
plied to the SFRs, since broad-band magnitudes are a 
poor tracer of the spatial extent of the [O II] line emis¬ 
sion. When we later discuss the relation between SFR 
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Input (Gyr) Input log 


Fig. 5. — Comparison between the input CSP maximum age and the output SSP age (left) or between the accumulated 

CSP mass and the SSP mass (right). The different symbols represent different values of the star formation timescale r (in 

units of Gyr) as reported in the left panel. The points in the left panel are given small random offsets to improve visibility. The gap at 
0-15 Gyr is caused by the sparse time sampling of our SSP models around this age. The diagonal lines represent the identity lines. 
The quantities and are largely consistent at >0.1 Gyr and r < a few 100 Myr, while becomes smaller at larger r due 

to the prolonged star formation. There is little systematic offset between the input and output masses. 


and stellar mass, we use the stellar mass within the fiber 
(^jfiber; stellar mass before the aperture correction) 
so that the two quantities refer to the same central part 
of the galaxies. Therefore, the present analysis does not 
map SFR and its relation to stellar mass to extended 
regions of the galaxies outside the fibers. 

3.3. Reliability of the present method 

Here the reliability of the present method is evaluated 
with Monte Carlo simulations. We generate 1,000 mock 
spectra by combining spectral models of quasar accretion 
disk, gas in the BLR, NLR, ISM, and host stellar pop¬ 
ulation as described in §3.1, with realistic PDFs for the 
parameter values. The power-law slope Opi is assumed 
to follow a Gaussian PDF with mean —1.5 and stan¬ 
dard deviation 0.3, which approximately reproduces past 
measurements (Ivezic et al. 2002; Pentericci et al. 2003). 
The SSP age is randomly drawn from 14 logarithmically 
spaced values between = 30 Myr and 10 Gyr, while 
the stellar mass is distributed uniformly on a logarithmic 
scale between log (M^/Mq) = 9 and 12. The ratio be¬ 
tween the mean amplitudes of the stellar model and the 
quasar power law is assumed to follow a logarithmically- 
uniform PDF between 0.01 and 10.0. The stellar velocity 
dispersion is drawn from a logarithmically-uniform PDF 
between 50 and 400 km s“^. All other model parameters 
are assumed to follow the Gaussian PDFs best describing 
the distributions derived in our own measurements. For 
each of the mock spectra, we randomly select a SDSS- 
RM quasar at z < 1 and apply its redshift and pixel flux 
errors to mimic a real spectrum. 

These mock spectra are processed through the de¬ 
composition algorithm used for the real data. Figure 
4 presents the comparisons between the intrinsic (in¬ 
put) and measured (output) parameter values. The stel¬ 
lar age and mass show relatively large scatter, reflect¬ 
ing the ambiguity in decomposing the nuclear and stel¬ 


lar components from the observed spectra. The scatter 
is significant at < 0.1 Gyr (the superscript “0” rep¬ 
resents the input parameter), where the relatively blue 
spectra of young stars are difficult to distinguish from 
those of quasar nuclei. These degeneracies result in over¬ 
estimated ages up to t* ~ 0.5 Gyr, hence we conserva¬ 
tively regard measured stellar ages at t* < 0.5 Gyr as 
upper limits (although we find only few such cases in the 
real quasars; see below). At larger ages, most of the scat¬ 
ter is caused by spectra with /* < 0.1 or A log M* > 
0.5 dex. When these simulated quasars are eliminated, 
as we did for the real sample, all six parameters are con¬ 
strained within a scatter of roughly 0.15 dex; the actual 
RMS scatter of each parameter is reported in each panel 
of the figure. Figure 4 also reports for each parameter 
the fraction (/id) of the simulated spectra whose out¬ 
put value is consistent with the input (i.e., the ordinate 
value is consistent with zero) within Icr statistical error. 
The fraction for cr* is low, due to the systematic over¬ 
estimation of this parameter at small cr° where the mea¬ 
surement is most difficult. For the other five parameters, 
/id are somewhat smaller than but roughly comparable 
to the expected 68 %. 

We also perform a set of simulations that replace 
the input SSP models with composite stellar population 
(GSP) models with exponentially declining star forma¬ 
tion history. CSP spectra are calculated as 

jCSP 

= f * - t)ff^it)dt; m = e-‘/^ 

^0 

(3) 

where t is the time since the onset of star formation 
(i.e., the maximum CSP age), ip(t) is the star formation 
history with the exponentially declining timescale r, and 
/f®^(t) is an SSP spectrum. We create 1,000 mock spec¬ 
tra by combining these CSP models with other emission 
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Fig. 6.— Comparison of M, (left; 140 successful measurements) or cr* (right; 73 successful measurements) for common objects measured 
by Shen et al. (2015b) and this work. The diagonal lines represent the identity lines. The two sets of measurements are in good agreement 
with each other, although there is a small (~10 km s“^) systematic offset in (j*. 


components as above, and process them through the de¬ 
composition algorithm (which uses the SSP models to 
represent stellar spectra). The results are shown in Fig¬ 
ure 5 for stellar age and mass; the other parameters are 
less sensitive to the choice of star formation history. The 
decomposition is unreliable for input ^0.1 Gyr, 

which is similar to the situation described above. At 
larger the output SSP age is consistent with 
when T is less than a few 100 Myr, while the former 
becomes systematically small at larger r due to the pro¬ 
longed star formation. There is little systematic offset 
between the output SSP mass and the input CSP mass, 
the latter being the mass accumulated during the whole 
star formation history. 

Finally, we compare our measurements of M, and cr* 
with the previous measurements by Shen et al. (2015b), 
who investigated the M, - cr* relation of the SDSS- 
RM quasars at 2 : < 1 using the same co-added spec¬ 
tra as in this work. They performed spectral decom¬ 
position with the principal component analysis (PCA) 
method following Vanden Berk et al. (2006). The PCA 
method has the advantage of being highly efficient in re¬ 
producing observed spectra with an optimal number of 
empirical eigenspectra, but the interpretation of the re¬ 
constructed spectra is usually not straightforward. The 
SMBH masses were derived with the Vestergaard & Pe¬ 
terson (2006) estimator as in this work, while the stel¬ 
lar velocity dispersion was measured with the vdispfit 
routine in the idlspec2d package and with the penal¬ 
ized pixel-fitting code (Cappellari & Emsellem 2004) in 
the spectral window A = 4125 - 5350 A. Figure 6 com¬ 
pares the M, and cr* values measured by the two studies. 
They are in good agreement with each other, although 
there is a slight systematic offset in cr*; our measure¬ 
ments give ~10 km s“^ larger values than those of Shen 
et al. (2015b) on average. The root-mean-square (RMS) 
scatter around this offset is ~ 40 km s“^, while the typ¬ 
ical measurement error is ~15 km s“^ in the both stud¬ 


ies.Our decomposition procedure is not optimized for 
the cr* measurements, since the spectral resolution of the 
adopted stellar library (cr ~ 80 km s“^) is not sufficiently 
high for this analysis; we chose this library since it covers 
the whole of our spectral window with a wide range of 
stellar parameters, which is crucial for a detailed study of 
the decomposed stellar emission. In contrast, the empha¬ 
sis is on the most reliable cr* measurements in Shen et al. 
(2015b), who used two dedicated routines adopting spec¬ 
tral libraries with high resolution (a < 25 km s“^). Our 
fitting procedure also includes Ca H and K lines, which 
tend to over-estimate cr* (Kormendy & Illingworth 1982; 
Bernardi et al. 2003; Greene & Ho 2006). 

4. RESULTS 

We now return to the actual spectra. The derived 
quasar and host properties are displayed in several pa¬ 
rameter planes in Figures 7 and 8. Their characteristic 
values (medians and 68% central intervals) are also pre¬ 
sented in Figure 9 (“Model 0”). The stellar ages cluster 
around t* = 1 Gyr and show a clear deficit at less than 
several hundred Myr, suggesting that the galaxies have 
not experienced major star formation in the recent past. 
The absence of young stars is consistent with the rela¬ 
tively low SFRs (within the spectroscopic fibers), which 
are less than 10 Mq yr“^ in most cases. Our results agree 
with those of Ho (2005) and Kim et al. (2006), who found 
that the SFRs in broad-line quasars at 2 : < 1 are ^ 10 
Mq yr“^ at most, and can be much smaller considering 
the contribution from AGN photoionization to lines like 
[O II]. On the other hand, Shi et al. (2009) report moder¬ 
ately higher SFRs (~ IOMq yr“^) than ours for a small 
number of SDSS quasars at 2 : ~ 1, based on observations 
of polycyclic aromatic hydrocarbon emission. At least a 
part of this discrepancy may be due to selection effects, 
as their sample is about an order of magnitude more lu¬ 
minous than the present sample. We stress that these 

Note that the two studies are not independent, since they use 
the same spectra for the cr* measurements. 
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are optically-selected broad-line quasars; dust-reddened 
or obscured quasars are different populations (see § 1 ) 
and may have much higher SFRs than those estimated 
here. The stellar and SMBH masses of the present sam¬ 
ple are in the range of IO^^Mq < M* < and 

IO^Mq < M, < 1 O®M 0 , respectively, suggesting that 
the SDSS-RM quasars are hosted by fairly massive galax¬ 
ies with massive SMBHs. Finally, Figure 7 demonstrates 
that objects with larger stellar fraction (/* > 0.5) tend 
to have lower T[oiii] and M, at a given redshift, while 
their host properties are not significantly different from 
those with smaller /*. 

We checked the robustness of our analysis to the model 
specifics, using the seven different basic models summa¬ 
rized in Table 1. Models 1 and 2 assume the metallicity 
Z = 0.5Zq and 2.0Zq, respectively, while Models 3 and 
4 adopt different IMFs from our default Chabrier (2003) 
IMF. In Models 5 and 6 , the host galaxies (i.e., the stel¬ 
lar and ISM model components) are subject to modest 
dust extinction, using the Small Magellanic Cloud type 
extinction curve (Rx = Ax/Eb-v = 1-39 (A/ 10000 )“^-^; 
Prevot et al. 1984; Richards et al. 2003; Hopkins et al. 
2004). We do not apply this additional dust extinc¬ 
tion to the quasar component, since the adopted empir¬ 
ical quasar templates already include the mean redden¬ 
ing, and broad-line quasars are usually observed with 
little extinction (see the discussion in §3.1). Finally, 
Model 7 assumes twice the broadening width of the Fe II 
template compared to the standard model. The spec¬ 
tral decomposition results with these different models 
are displayed in Figure 9. The lower/higher metallic¬ 
ity leads to older/younger stellar ages, as expected from 
the well-known age-metallicity degeneracy, and slightly 
higher/lower stellar masses. The Salpeter (1955) IMF 
results in ~0.2 dex higher M* than the Chabrier (2003) 
or Kroupa (2001) IMF, which reflects higher mass-to- 
luminosity ratios by a similar factor. Assuming a color 
excess of Eb-v = 0.2 or 0.5 mag results in the stellar 
masses and SFRs simply increasing by the correspond¬ 
ing extinction factors {Ay ~ 0.6 or 1.5 mag), while 
the other parameters remain nearly unchanged. The 
assumption of the broadening width of the Fe II tem¬ 
plate has little impact. Overall, our results are qualita¬ 
tively insensitive to the specifics of the spectral models; 
the quasar host galaxies are characterized by intermedi¬ 
ate stellar ages (t* ~ 1 Gyr), fairly large stellar masses 
(M* ^ IO^^Mq), relatively low SFRs (< IQMq yr“^ un¬ 
less the dust extinction is very large), and large SMBH 
masses (M, ~ lO^M©). 

As an additional test of the reliability of the SFR es¬ 
timates, we compare them with SFRs estimated with a 
different calibrator, using the rest-frame M-band luminos¬ 
ity (L„) of the stellar component. We use the Hopkins 
et al. (2003) calibration to derive the Lu-based SFR: 


SFR (Mq yr-i) 




1.186 


1.81 X 1021 W Hz-1 


(4) 


with Lu calculated from the decomposed host spectrum 
of each object, extrapolated to A < 3700 A (i.e., outside 
the fitting window) with the best-fit stellar model. The 
[O II]-based and T„-based SFR estimates are subject to 
different uncertainties. As described in §3.2, the [O II]- 


TABLE 1 
Model specifics. 


Model 

Description 

U“ ... 

standard 

1 ... 

Z = 0.5Zq 

2^ ... 

Z = 2.0Zq 

3 ... 

Kroupa (2001) IMF 

4 ... 

Salpeter (1955) IMF 

5 ... 

Eb-v = 0.2 mag 

6 ... 

Eb-v = 0.5 mag 

7 ... 

crpell = 3,000 km 


®The standard model adopts Z = Zq, Chabrier (2003) IMF, 
Eb-v = 0.0 mag, and (TFell = 1,500 km s”!. 

I’The STELIB stellar library is replaced with MILES which has 
a wider coverage of stellar ages at this metallicity in the Maraston 
& Stromback (2011) models. 

based SFRs may be affected by object-to-object variation 
of the quasar contribution. The dust extinction affects 
nebular lines and stellar continuum in different ways, due 
to the different geometry of the emission regions (Calzetti 
et al. 1994; Chariot & Fall 2000). Nonetheless, Figure 10 
shows that the two estimates are roughly comparable, al¬ 
though there is a large scatter. It may not be surprising 
that dust has only limited effect on these estimates. As 
previously mentioned, past observations demonstrated 
that the lines of sight toward optically-selected broad¬ 
line quasar nuclei are almost dust free; the color excess 
is less than Eb-v = 0-1 mag in the vast majority of 
cases (Richards et al. 2003; Hopkins et al. 2004; Salvato 
et al. 2009; Matute et al. 2012; Krawczyk et al. 2014). 
A galaxy-wide dust screen is not likely to be present, 
although there could be localized dusty star-forming re¬ 
gions. However, the present results indicate that the host 
galaxies are dominated by stars with intermediate ages 
and are not active star formers (at least within the spec¬ 
troscopic fibers). Specific SFR (SFR/M,) is even less 
sensitive to dust extinction than is SFR. Our test indi¬ 
cates that assuming moderate extinction in the decom¬ 
position procedure (Models 5 and 6 ) results in the stel¬ 
lar mass estimates simply increasing by the correspond¬ 
ing extinction factors. This is because the extinction 
does not alter the derived stellar age and hence mass- 
to-luminosity ratio significantly. As a result, extinction 
factors in SFR and M* are effectively cancelled out when 
specihc SFR is derived. 

We emphasize that our sample is not volume limited. 
However, it includes most of the SDSS-RM quasars at 
z < 1, which are expected to be fairly complete down to 
the limiting magnitude (see §2). Most of the initial sam¬ 
ple is successfully processed through the decomposition 
analysis; the fraction of objects whose results we deem 
unreliable (/* < 0.1 or A log M* >0.5 dex) is only 18 
%. We provide the decomposition results for all the 191 
quasars in the electronic edition of the journal. 

5. DISCUSSION 

5.1. Host galaxies of unobscured quasars 


Little effect of dust extinction on the derived stellar age is ex- 
pected from the narrow spectral window used in the first step of the 
decomposition procedure (A = 3700 — 4050 A; see §3.1), in which 
the extinction cannot alter the spectral slope significantly. How¬ 
ever, the present argument may be complicated by the fact that 
nebular lines and stellar continuum are affected by the extinction 
in different ways, as described above. 
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Fig. 7. — Spectral decomposition results of our sample of 156 quasars with successful measurements. Shown are (top left), M* (top 
right), SFR (middle left), cr* (middle right), F[oiii] (bottom left), and M, (bottom right), all as a function of redshift. The objects with 
the stellar fraction /* > 0.5 are marked with the large open circles. The error bars indicate statistical errors; there may be additional 
systematic uncertainties as discussed in the text. The dotted line in the top left panel represents the minimum age of our SSP models, t* 
= 30 Myr; the two objects on this line are well fitted with this age but may have even younger stellar populations. 


Before discussing the results, we first compare our 
measurements with those of Matsuoka et al. (2014a). 
They studied host galaxies of SDSS broad-line quasars 
at z < 0.6 with the imaging decomposition technique 
(see §1), which uses completely different data and meth¬ 
ods from those in the present analysis. The two studies 
share no common objects. We calculated the host col¬ 
ors of the present sample by convolving the decomposed 
galaxy spectra (in the rest frame; they are extrapolated 
to A < 3700 A, i.e., outside our fitting window, with the 
best-fit stellar models) with the SDSS filter transmissions 
and created the CMD as shown in Figure 11. As a ref¬ 
erence, we also plot the distribution of inactive galaxies 
at 0.5 < z < 1.0 taken from the COSMOS/UltraVISTA 
AT-band selected catalog (Muzzin et al. 2013); the rest- 
frame u — r colors were calculated with the best-fit spec¬ 
tral model of each object given in the catalog. The two 


measurements of the quasar hosts broadly agree in distri¬ 
butions of stellar mass and color, demonstrating that the 
quasars are preferentially hosted by massive galaxies dis¬ 
tributed from the massive tip of the blue cloud to the red 
sequence. The median stellar age of the present sample, 
t* 1 Gyr, corresponds to the gap of the bimodal color 
distribution of inactive galaxies. It is intriguing that the 
quasars occupy the colors centered on this gap, where 
star-forming galaxies may be rapidly transitioning to the 
quiescent phase due to the quenching of star formation. 
The stellar masses of the quasar hosts (M* ~ IO^^Mq) 
broadly agree with the turn-over mass of the local galaxy 
mass function at the high mass end (e.g., Bell et al. 
2003). Galaxy growth is significantly suppressed above 
this mass (Peng et al. 2010), for which quasar activity 
may be playing an important role. 

Figure 12 displays the relation between SFRs and stel- 
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log(Mstar/Msun) log(MBH/M3,J log (L[o„i]/erg s"’) 

Fig. 8.— Relations between (top panels) or SFR (bottom panels) and M* (left), M* (middle), or Ffoiii] (right). All the successful 
measurements of the 156 quasars are plotted. The objects with the stellar fraction /* > 0.5 are marked with the large open circles. The 
dotted lines in the top panels represent the minimum age of our SSP models, f* = 30 Myr. 


lar masses within the spectroscopic fibers/"' along with 
the so-called “main sequence” of star-forming galaxies 
(e.g., Brinchmann et al. 2004; Noeske et al. 2007; Daddi 
et al. 2007; Elbaz et al. 2007, 2011; Steinhardt et al. 
2014). The ratio between these two quantities defines 
specific SFR, which is not very sensitive to the dust ex¬ 
tinction in the present analysis as discussed in §4. The 
majority of the quasar hosts fall below the main sequence 
at 2 ; = 1, i.e., their star formation efficiencies are consid¬ 
erably lower than those of normal star-forming galax¬ 
ies at similar redshifts. On the other hand, SFRs mea¬ 
sured in host galaxies of dust-obscured. X-ray selected 
AGNs at similar redshifts tend to be higher than esti¬ 
mated here and trace the main sequence of star-forming 
galaxies (Mullaney et al. 2012; Santini et al. 2012). How¬ 
ever, Mullaney et al. (2015) demonstrated recently that 
these SFR estimates are biased toward higher values due 
to bright outliers in stacking analysis of infrared data, 
and that the actual SFR distribution of obscured AGNs 
is broader than, and has ^0.4 dex lower peak of, that of 
main-sequence galaxies. These findings are broadly con¬ 
sistent with our results for broad-line quasars presented 
in Figure 12. 

The intermediate stellar ages and relatively low SFRs 
of the present sample suggest that the quasar hosts are 
linked to post-starburst populations. Since the first dis- 

Here we use the stellar masses within the fibers (i.e., before 
the aperture correction) so that both stellar masses and SFRs refer 
to the same central region of the galaxies; see §3.2. 


covery by Dressier & Gunn (1983), post-starburst galax¬ 
ies have been recognized as an important stage of galaxy 
evolution. They are characterized by the presence of 
strong B aimer absorption lines and the absence of emis¬ 
sion lines, which are commonly interpreted as evidence 
for starburst activity within the past ~1 Gyr, which was 
quenched abruptly (e.g., Quintero et al. 2004; Goto 
2005). The cause of this quenching is still unknown, 
but the present study indicates that the formation of at 
least some of post-starburst galaxies is related to the nu¬ 
clear activity. Quasars detected in post-starburst galax¬ 
ies (e.g., Brotherton et al. 1999, 2002) may represent the 
relevant phase. In the Hopkins et al. (2006) merger- 
driven co-evolution model, a galaxy experiences vari¬ 
ous transition phases after a gas-rich merger, including 
a luminous infrared galaxy phase with intense star for¬ 
mation and dust production, an obscured AGN/quasar 
phase with activated central SMBHs, and an unobscured 
AGN/quasar phase before settling into a quiescent early- 
type galaxy. Our results point to the link between 
quasar hosts and post-starburst galaxies, suggesting that 
merger-induced starbursts are largely quenched before 
the dust obscuring the nuclear region is cleared away 
and the central quasar becomes optically visible. Re¬ 
cently Pattarakijwanich et al. (2014) estimated the num¬ 
ber density of post-starburst galaxies to be roughly 10“® 
Mpc“'^ mag“^ for M 5000 > —20 mag at 0 < z < 1, where 
Afsooo is the absolute magnitude at rest-frame 5000 A. 
This value is roughly comparable to the number density 
of broad-line quasars at similar redshifts and luminosities 
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Model Model 

Fig. 9.— Spectral decomposition results with different model 
specifics as summarized in Table 1; Model 0 (standard), 1 (Z = 
O. 5 Z 0 ), 2 {Z = 2.0Zq), 3 (Kroupa IMF), 4 (Salpeter IMF), 5 
{Eb-v = 0-2 mag), 6 (Eb-v = 0.5 mag), and 7 (fTpeii = 3000 km 
s~^). The median values (circles) and 68 % central intervals (error 
bars) of the derived parameter distributions are shown for (top 
left), M* (top right), SFR (middle left), cr+ (middle right), T[oiii] 
(bottom left), and M* (bottom right). 



[Oll]-based SFR yr ') 

Fig. 10.— Comparison between the L^-based and [O Il]-based 
SFRs. The diagonal solid line represents the identity line. The 
two estimates are roughly comparable, although there is a large 
scatter. 


(Richards et al. 2006; Groom et al. 2009). Kauffmann et 
al. (2003) noted that the strong absorption observed 
in low-z {z ~ 0.05) SDSS narrow-line AGNs indicates 
they experienced a burst of star formation that ended 
within the past ~1 Gyr. Yesuf et al. (2014) demonstrated 
that the AGN fraction is three times higher in post- 
starburst galaxies than in inactive galaxies. However, 
they also reported that (obscured) AGN phase is delayed 


by a few 100 Myr relative to post-starburst phase, and 
hence cannot play a primary role in the quenching of 
starbursts. 

There have been a number of studies comparing the 
color and mass distributions of AGN host galaxies with 
those of inactive galaxies (e.g., Nandra et al. 2007; Salim 
et al. 2007; Hickox et al. 2009; Schawinski et al. 2010). 
Recent work suggests that, using proper mass-limited 
samples, the AGN fraction is constant as a function 
of host color, or may be enhanced in blue star-forming 
galaxies (e.g., Silverman et al. 2009; Aird et al. 2012; 
Hainline et al. 2012; Rosario et al. 2013). Luminous 
AGNs are preferentially found in massive host galax¬ 
ies, which is perhaps due to the presence of massive 
SMBHs at the centers and/or of large reservoirs of gas 
in such galaxies. Most of these studies are based on ob¬ 
scured AGNs selected via their X-ray luminosity or opti¬ 
cal narrow emission lines. Our work (Figure 11) suggests 
that broad-line quasar hosts prevail in the color range 
1.5 < M — r < 2.5 above masses of log (M^/Mq) ~ 10.5 
and that the quasar fraction increases toward bluer col¬ 
ors, a trend that continues to the end of the galaxy color 
distribution at u — r < 1.5. In order to show this trend 
more clearly, we create a mass-matched sample of inac¬ 
tive galaxies; for each of the quasar hosts, the galaxy 
with the closest mass and redshift are drawn from the 
GOSMOS/UltraVISTA AT-band selected catalog. Figure 
13 compares the histograms of the rest-frame (u — r) col¬ 
ors of the two samples. The quasar host distribution is 
shifted toward bluer colors than that of mass-matched 
inactive galaxies, demonstrating that bluer galaxies are 
more likely to host quasars. 

Meanwhile, a drawback of the present analysis is the 
lack of far-infrared information and hence the dust- 
absorbed part of star formation. Although a galaxy-wide 
dust screen is unlikely to be present (see §4), the present 
work cannot trace localized star-forming regions hidden 
by dust at UV-to-optical wavelengths. Such regions can 
be mapped reliably only with far-infrared observations, 
which we will exploit in future work. Our analysis is also 
limited to within the spectroscopic fibers, whose 2" di¬ 
ameter corresponds to approximately 14 kpc at the me¬ 
dian redshift of the sample (see §3.2). If the majority 
of the star formation took place at outside this aperture, 
then our results would become inconsistent with the past 
studies which measured the entire host galaxies. 


5.2. Scaling relations 

Our analysis provides both stellar velocity dispersions 
and masses as well as SMBH masses, which allows us to 
investigate the scaling relations between them. Figure 
14 presents the M, — and M, — M* relations of the 
present sample. We calculate the linear regression lines 
with the IDL routine MPFITEXY, as: 


log 


M. \ 
Mq) 


and 


log 


M. 

IQSMq 


= (1.64±0.29) logf- — -r ) 

^ ^ ^V200 kms”V 

±(0.43 ±0.05), 

) = (0.36 ± 0.08) log (^) 

±(0.21 ±0.04). 


(5) 
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Fig. 11.— Rest-frame {u — r) colors and stellar masses of the quasar host galaxies in this work (red dots) and in Matsuoka et al. 
(2014a, blue dots). The typical error for the present sample is shown by the error bar at the bottom right corner. The small gray dots 
and contours, drawn at logarithmically stepped levels of number density, represent non-AGN galaxies at 0.5 < z < 1.0 taken from the 
COSMOS/UltraVISTA A'-band selected catalog (Muzzin et al. 2013). The dots are given small random offsets to improve visibility. The 
dashed lines mark the SSP colors with t* = 0.03, 0.1, 0.3, 1, 3, and 10 Gyr, as labelled. The quasar hosts are preferentially hosted by 
massive galaxies distributed from the massive tip of the blue cloud to the red sequence. Their mean stellar age (t* ~ 1 Gyr) corresponds to 
the gap of the bimodal distribution of inactive galaxies, where blue star-forming galaxies may be rapidly transitioning to the red sequence. 


The intrinsic scatters are 0.39 and 0.44 dex for the 
M,—cr* and M, — M* relations, respectively. The M,—ct* 
relation is almost identical to that of Shen et al. (2015b, 
see §3.3), which is not surprising given the consistent 
measurements of the two quantities (Figure 6). Our sam¬ 
ple shows no evolution of the mean M./cr* ratio from the 
local Universe, while the slope of the relation appears to 
flatten toward high redshifts. This flattening was also 
reported by Shen et al. (2015b) and was interpreted as a 
result of a selection bias; with simple simulations, they 
concluded that a non-evolving M, — ct* relation is favored 
at z < 1. 

We also found similar trends in the M, — M* relation. 
Our Mif estimates include a galaxy disk component if it 
is present, hence they should be regarded as the upper 
limits of the bulge masses (M*^buige)- A similar compar¬ 
ison between the M, — M* relation at z ~ 1.4 and the 
local M, — M*_buige relation was presented by Jahnke et 
al. (2009), who found no significant difference between 
the two. Since their sample shows evidence for substan¬ 
tial disk components, the authors concluded that bulges 
have grown predominantly via redistribution of the disk 
mass since z = 1.4 to the present epoch. This conclu¬ 


sion was supported by Sun et al. (2015) with a larger 
sample at similar redshifts and with a careful assessment 
of selection biases. Our results are consistent with these 
findings, but we have no estimate of disk contribution to 
M^, and cannot provide a direct support for the above 
scenario. The present M, — M* relation is also similar to 
that of Matsuoka et al. (2014a), who reported that the 
relation is offset toward large M, compared to the local 
relation at fixed M*. This result is likely due to the shal¬ 
low slope of the relation, producing systematically large 
M,/M* ratios at M* < IO^^Mq. 

5.3. Applicability of the present analysis 

The present spectral decomposition analysis is made 
possible due to the low-to-moderate luminosities of the 
quasars, with relative stellar contribution /* > 0.1 mea¬ 
sured at A = 4000 A, and the high SNR of the spec¬ 
tra. Meeting both conditions usually requires significant 
telescope resources; in this work, we have exploited the 
co-added SDSS-RM spectra whose net exposure times 
amount to 65 hours per object. As we demonstrated in 
the previous sections, clear detection of stellar absorption 
features allows for reliable extraction of the host com- 
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Fig. 12.— Relation between SFRs and stellar masses within the 
2”-diameter spectroscopic fibers. All the successful measurements 
of the 156 quasars are plotted. The solid and dotted lines represent 
the main sequence of star-forming galaxies and their 68% confi¬ 
dence intervals at 2 = 1 and 2 = 0, respectively, taken from Elbaz 
et al. (2007). The quasar hosts fall below the main sequence, i.e., 
their star formation efficiencies are considerably lower than those 
of normal star-forming galaxies at similar redshifts. 



Fig. 13. — Histograms of the rest-frame (u—r) colors of the quasar 
host galaxies in this work (thick line) and the mass-matched inac¬ 
tive (non-AGN) galaxies drawn from the COSMOS/UltraVISTA 
catalog (thin line). The quasar host distribution is shifted toward 
bluer colors than that of inactive galaxies, demonstrating that bluer 
galaxies are more likely to host quasars. 

ponents from the observed spectra. In this section, we 
give a rough guideline for estimating whether the present 
analysis is applicable to a given spectrum. 

The first requirement is a certain fraction (/* > 0.1) 
of host galaxy light in the total spectrum. This fraction 
can be estimated roughly from the depths of the stellar 
absorption lines. Figure 15 presents the stellar fraction 
/* as a function of EW (Ca K), the rest-frame equiv¬ 
alent width of the Ca K absorption line. Ca K is one 
of the strongest stellar absorption features in a quasar 
spectrum at UV-to-optical wavelengths, while Ca H is 


less visible due to the superposed [Ne III] A3969 emis¬ 
sion. We measure EW (Ca K) in the wavelength range A 
= 3920 - 3955 A, with the blue and red continuum levels 
defined at A = 3880 - 3920 A and A = 3990 - 4030 A, 
respectively. Figure 15 indicates that EW (Ca K) and /* 
are moderately correlated. The majority of the quasars 
with EW (Ca K) > 2 A have more than 10% of the light 
at 4000 A from the host, and stellar emission is the major 
contributor to a quasar spectrum on average when one 
observes EW (Ca K) > 6 A. 

The second requirement is high SNR. We test the min¬ 
imum SNR required to perform a reliable decomposition 
analysis by artificially degrading the SDSS-RM data and 
determining how the derived uncertainties change as the 
SNR decreases. We take the original high-SNR spectra 
with successful decomposition and add noise assuming 
a Gaussian PDF with standard deviations of 2, 4, 8, or 
16 times the original pixel errors. These mock spectra 
are processed through the decomposition algorithm. We 
omit the calibration errors in Equations 1 and 2 in these 
simulations to view the effect of decreasing SNR directly. 
Figure 16 displays the resultant uncertainties in the six 
physical quantities as a function of the SNR per pixel 
measured at A = 4000 A. These plots can be used to 
infer the required SNR to achieve the desired accuracy 
in one of the parameters; for example, a rough estimate 
of stellar mass (Alog M* < 0.5 dex) is possible at SNR 
3 on average, while stellar age is poorly constrained 
(At,/t» ^ 1) at that SNR. All six parameters are mod¬ 
erately well constrained at SNR > 10, where the actual 
SDSS-RM quasars are located. 

6. SUMMARY 

We present the results of a spectral decomposition 
analysis of broad-line quasars at z < 1. The sample is 
drawn from the SDSS-RM project, in which a single spec¬ 
troscopic field was repeatedly observed with the BOSS 
spectrograph to explore the variability of quasars. The 
high SNR co-added spectra of 191 SDSS-RM quasars, 
whose net exposure times amount to approximately 65 
hours per object, are used for the analysis. We develop a 
technique to decompose a quasar spectrum into nuclear 
and host components using spectral models of the quasar 
accretion disk, gas in the BLR, NLR, ISM, and stars in 
the host galaxy. This procedure provides estimates of 
stellar age and mass, SFR, stellar velocity dispersion, 
quasar [O III] luminosity, and SMBH mass, among other 
quantities. The decomposition was successful for more 
than 80 % of the initial sample; these objects have more 
than 10 % of their light at 4000 A from the host, and the 
error in stellar mass estimates Alog (M^/Mq) is less than 
0.5. We find that the quasars are preferentially hosted 
by massive (M* ~ IO^^Mq) galaxies with characteris¬ 
tic stellar ages around t* = 1 Gyr. It coincides with the 
gap of the bimodal color distribution of inactive galaxies, 
where blue star-forming galaxies may be rapidly transi¬ 
tioning into the red sequence. The deficit of young stars 
is consistent with the estimated low SFRs, which place 
the quasar hosts below the main sequence of star-forming 
galaxies at similar redshifts. These facts suggest that the 
quasar hosts have experienced an episode of major star 
formation sometime in the past ^1 Gyr, which was sub¬ 
sequently quenched or suppressed (at a slower rate) . We 
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Fig. 14.— Relation between M* and (j* (left) and between M# and Af* (right). The thick solid lines represent the regression lines of the 
present sample, while the thin solid lines indicate those of the previous measurements by Shen et al. (2015b, “S15”) and Matsuoka et al. 
(2014a,b, “M14”). The dotted and dashed lines represent the local relations taken from Kormendy & Ho (2013, “KH13”) and Haring &; 
Rix (2004, “HR04”), respectively. Note that the KH13 and HR04 relations have been calibrated with bulge masses, while M* of this work 
includes a galaxy disk component if it is present (see text). The mean M#/cr* or M*/M* ratio of the quasar hosts is consistent with the 
local value, while the slope of the relation appears to flatten toward high redshift. This flattening can be explained by selection biases (see 
text). 



Fig. 15. — Stellar fraction (/*) measured at A = 4000 A in our 
sample of 191 quasars, as a function of the rest-frame EW (Ca K); 
a positive EW value means absorption. The large circles represent 
median values in bins of EW (Ca K). The dotted line marks the 
success threshold of our decomposition analysis, /♦ > 0.1. There 
is a moderate correlation between the two quantities, which allows 
one to make a rough estimate of /* by measuring EW (Ca K) in a 
given spectrum. 

also find that the scaling relations between M, and cr* 
and between M, and M* are consistent with our previ¬ 
ous measurements. The mean or M,/M* ratio is 

consistent with the local value, while the slope of the rela¬ 
tion is flatter than the local relation. As demonstrated by 
Shen et al. (2015b), this result is likely due to a selection 
bias, and no evolution of the scaling relations is favored 
at z < 1. Finally, we demonstrate that the stellar contri¬ 
bution to a quasar spectrum can be inferred roughly from 
the equivalent widths of Ca K absorption line, and pro¬ 


vide estimates of the minimum SNR required to perform 
a spectral decomposition with desired precision. 

The present work will be greatly advanced with future 
powerful spectrographs such as the Subaru Prime Fo¬ 
cus Spectrograph (PFS; Takada et al. 2014), which will 
produce large numbers of high-quality spectra of galax¬ 
ies and quasars out to high redshifts. Our method, op¬ 
timized for BOSS fiber spectra, can be readily applied 
to other high-quality spectra in order to study quasars 
with lower luminosity and higher redshifts. It is also im¬ 
portant to extend the present analysis into longer rest- 
frame wavelengths (A > 5400 A) with near-infrared spec¬ 
trographs, where the relative host contribution becomes 
larger and the possible effects of dust extinction become 
smaller. 
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Fig. 16.— Uncertainties in t* (top left), M* (top right), SFR, (middle left), a* (middle right), I/[oiii] (bottom left), and M% (bottom 
right), as a function of the SNR per pixel measured at A = 4000 A. The large dots represent the original SDSS-RM spectra, while the 
small dots represent the mock spectra with noise added. The quantity has a lower limit at 0.2 (dotted line in the upper left panel) 

because of the 20 % error added to account for non-contiguous age coverage of the SSP models (see §3.2). The squares connected by the 
solid lines show the median values in bins of SNR. 
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